Quasiparticle decay rate of Josephson charge qubit oscillations 
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We analyze the decay of quantum oscillations in a charge qubit consisting of a Cooper pair 
box connected by a Josephson junction to a finite-size superconductor. We concentrate on the 
contribution of quasiparticles in the superconductors to the decay rate. Passing of a quasiparticle 
through the Josephson junction leads to the escape of the qubit out of its Hilbert space, and thus 
determines the decay rate of quantum oscillations. We find the temperature dependence of the 
quasiparticle contribution to the decay rate for open and isolated systems. The former case is 
realized if a normal-state trap is included in the circuit, or if just one vortex resides in the qubit; 
we find exponential suppression of the rate, F oc exp (— A/T), at low temperatures (here, A is the 
superconducting gap). In a superconducting qubit isolated from the environment F oc exp (— 2A/T) 
if the number of electrons is even, while for an odd number of electrons the decay rate remains finite 
in the limit T — » 0. We estimate F for realistic parameters of a qubit. 

PACS numbers: 03.65.Yz, 03.67.Lx, 74.50.+r, 85.25.Cp 



I. INTRODUCTION 

Recent experiments with superconducting qubits show 
promising results, allowing one to observe hundreds of 
qubit charge or phase oscillation»i*2iiii. A superconduct- 
ing qubit is essentially a controllable quantum two-level 
system and can be realized using phase or charge degrees 
of freedom. Most of the charge qubits use a Cooper- 
pair box (CPB), a small mesoscopic island connected 
to a large superconducting reservoir via two Josephson 
junctions. A device with a large superconducting gap, 
A > E c > Ej ^> T, can be controlled with the gate volt- 
age and magnetic flux, and has only one discrete degree 
of freedom: number of Cooper pairs in the box. (Here, 
E c is the charging energy of the island, and Ej is an ef- 
fective Josephson energy of its junctions with the reser- 
voir; Ej depends on the flux.) At the degeneracy point 
the qubit evolves coherently between the ground and ex- 
cited states, |— ) and |+), which are represented by the 
superposition of or 1 excess Cooper pairs in the box. 
The frequency of these oscillations is determined by the 
Josephson energy. 

The main difficulty in technological realization of a 
charge-based superconducting qubit is due to the deco- 
herence present in it. The mechanisms of the decoher- 
ence are currently unknown and can be attributed to 
phonons, two-level systems in insulating barrier, ther- 
mally excited quasiparticles*^. In this paper we con- 
centrate on the contribution of the quasiparticles to the 
decay rate. Quasiparticles having a continuous spec- 
trum are inherently present in any superconducting de- 
vice and set a fundamental constraint on the coherence 
time. Quasiparticle "poisoning" , first investigated in the 
context of charge-parity effects in mesoscopic supercon- 
ductorsi^, manifests itself also in the experiments with 
charge qubitsii*i2ii&i2ii£. It was reported that even at 
low temperatures (~10-50mK) quasiparticles are present 
in the CPB. If this is the case, Hilbert space of the CPB 



expands, and the qubit is no longer a simple two-level 
system. The transient presence of a quasiparticle in the 
CPB detunes the qubit from the resonant state of Cooper 
pair tunneling and affects coherent oscillations. In this 
article we build a quantitative theory of the quasiparticle 
effect on the charge qubit oscillations. 

We consider two regimes which can be realized exper- 
imentally: open system corresponding to fixed chemi- 
cal potential in the reservoir, and isolated system corre- 
sponding to fixed number of electrons in the qubit (see 
Fig.l). The former case allows for a change of the to- 
tal number of electrons in superconducting parts of the 
system, and is experimentally realized if a quasiparticle 
trap, e.g. normal-state part or a vortex, is included in the 
circuit. The latter case corresponds to a superconducting 
qubit isolated from the normal-metal environment. Both 
cases may be relevant in the context of the cavity quan- 
tum electrodynamics experiments where the state of the 
qubit is determined using photon degrees of freedom. 

The paper is organized as follows. We begin in Sec. 
II with a brief overview of the charge qubit before the 
discussion of the quasiparticle effect on the charge qubit 
in an open system. In Sec. Ill we consider the opposite 
case of the isolated qubit. Finally, in Sec. IV we present 
simplified results for the quasiparticle contribution to the 
decoherence rate for the mentioned experimental realiza- 
tions and discuss how to decrease it. 



II. STATES OF THE QUBIT IN AN OPEN 
SYSTEM 



Dynamics of the superconducting charge qubit is usu- 
ally described by an effective Hamiltonian 



H eff =E c (N-N g ) 2 + H J , 



(1) 



where E c is charging energy of the box, TV is di- 
mensionless charge of the CPB, and N g is dimension- 
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A. Thermodynamic properties of a qubit in an 
open system 

At a finite temperature, density of quasiparticles is ex- 
ponentially small, n oc exp (— A/T), but the number of 
quasiparticles in the grain can be of the order of one. It 
is important to point out that even one unpaired electron 
can affect the qubit performance. In order to estimate 
the number of quasiparticles in the grains, one has to 
account for the huge statistical weight of the states with 
a single quasiparticle, proportional to the volume of the 
grain, N ~ V '2-k AbTvbVb exp (— A&/T), where Vb is nor- 
mal density of states(per volume), Vb volume of the grain. 
For an isolated dot, the characteristic temperature^ at 
which quasiparticles appear is (fee = 1) 



n = 



As 



H^b/Sb) ■ 



(3) 



FIG. 1: Schematic picture of the charge qubit in different ex- 
perimental regimes: a) open system, b) isolated system. The 
left superconducting mesoscopic island is the Cooper pair box 
connected via a tunable Josephson junction to the large su- 
perconducting reservoir (right). Gate bias is applied through 
the capacitance C a (assuming that Coo C 9 ). 



less gate voltage. The Hamiltonian corresponding to 
tunneling of Cooper pairs, Hj, is defined as Hj = 
(\N + 2) (N\ + h.c), where Ej is effective Joseph- 
son energy. In the regime where superconducting gap is 
the largest energy scale in the system A > E c > Ej ^> T, 
quasiparticles are usually neglected, and the dynamics 
of the system is described by the above Hamiltonian, 
where there is only one discrete degree of freedom - ex- 
cess number of Cooper pairs in the box. At the operat- 
ing point, when the dimensionless gate voltage is tuned 
to be equal to one, N a = 1, there is degeneracy with re- 
spect to charging energy between the charge states \N) 
and |JV + 2). This degeneracy is lifted by the Joseph- 
son energy, and the states of a qubit are described by 
the symmetric and antisymmetric superposition of the 
charge states, R = |Af+ ™ and |+) = ^ 
with corresponding energies: 



where 8b 



is the mean level spacing in the box, 



E c and uj + = E c 



El 

2 ' 



(2) 



Other charge states have much higher energy, and effec- 
tively the Cooper pair box reduces to a two-level sys- 
tem. Once a qubit is excited, quantum oscillations be- 
tween states |— ) and |+) emerge, and the frequency of 
these oscillations is determined by the Josephson energy 
EjlQ. The appearance of a quasiparticle with a contin- 
uum spectrum provides a channel for relaxation of the 
qubit. Since quasiparticles are inherently present in any 
superconducting system, their contribution to the decay 
rate is intrinsic. 



and Ab is superconducting gap energy in the box. The 
appearance of a quasiparticle in the qubit occurs at lower 
temperature T b * due to the finite charging energy: 



2? 



A b -E c + Ej/2 
\og{A b /Sb) 



t: i 



E c - E.j/2 



(4) 



If T <C T b * , states with odd number of electrons in the 
box are statistically rare. The probability of finding the 
qubit in a "good" state (not poisoned by quasiparticles) 
is important for qubit preparation and is determined by 
thermodynamics. However, the qubit coherence time is 
controlled by kinetics, which we study in the next section. 



B. Quasiparticle decay rate in the open system 

If an unpaired electron tunnels into the CPB, it tunes 
the qubit away from the resonant state of Cooper pair 
tunneling, which leads to the decay of quantum oscil- 
lations. At the operating point N g = 1 it is energeti- 
cally favorable for a quasiparticle to tunnel to the CPB, 
since charging energy is gained in such process. Assum- 
ing that initially the qubit was prepared in the state with 
no quasiparticles in the box, the lifetime of the qubit is 
determined by the time of quasiparticle tunneling to the 
CPB. In order to estimate this time we use the following 
Hamiltonian: 



H = H n 



— Hbcs 



H 



BCS 



H 7 



E c (Q/e - N g ) 2 



(5) 



where H% cs , 



Hg CS are BCS Hamiltonians of the box 



and superconducting reservoir (see Fig. 1), and Q de- 
notes the charge in the box. The tunneling Hamiltonian 
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between two electrodes, Ht, is denned as 



,t 



H.c.) 



(6) 



kpa 



where ife p is tunnelling matrix clement, c kM , c PiCr are the 
annihilation operators for an electron in the state k,a in 
the CPB and state p, a in the superconducting reservoir. 

We now consider the lifetime of the qubit in the open 
system, allowing for a change in the number of electrons 
in the superconducting reservoir. Assuming that the 
qubit was prepared in the initial state without a quasi- 
particle in the box, the time of its coherent evolution 
is limited by the rate of quasiparticle tunneling to the 
CPB. In order to calculate the lifetime of the qubit, we 
have to distinguish between tunneling of Cooper pairs 
and quasiparticles. To do this, we write the Hamiltonian 
H in Eq. JSJ in the following form: 



H = Hn 



Hi, 



(7) 



where Hq = H + Hj, and Hi = Ht — Hj. Hj is second 
order in tunneling amplitude 

1 



Hj = \N) (N\ H T 



E-H, 



■H T \N+2) (N+2\ + H.c. (8) 



The matrix element (N\ Httt-7F H t \ n + 2 > is propor 

El -On 



tional to effective Josephson energy Ej, and Ht is de- 
fined in Eq. JfJJl. Without quasiparticles Hamiltonian Hq 
reduces to Eq. (JTJ. 

The quasiparticle tunneling rate is found using Fermi's 
golden rule and averaging over initial configuration with 
the appropriate density matrix (% = 1) 



2tt 



{flHiMfSiEf-EJptfHo) 



(9) 



Here, p(/3Hg) is the density matrix for the initial state 
of the system. The perturbation Hamiltonian accounts 
for quasiparticle tunneling only, (f\Hi\i) = (f\H T \i). 
Thus, Eq. © takes into account Cooper pair tunneling 
exactly while treating quasiparticle tunneling perturba- 
tively. At the operating point, N g = 1, initial state of 
qubit is defined by the superposition of and 1 excess 
Cooper pairs in the box, |±) = l^liL^±2lj the final state 

|/) = |7V + 1) is the state with odd number of electrons 
in the CPB, corresponding to charge le. There are two 
mechanisms that contribute to the rate of the process 
|±) — * \N + 1): (1) a quasiparticle tunnels from the su- 
perconducting reservoir to the CPB, and (2) a Cooper 
pair in the box breaks into two quasiparticles, and then 
one quasiparticle tunnels out into the reservoir. Two cor- 
responding contributions to the total tunneling rate are 

r ± = r 1± + r 2± , (10) 

T 1± = 2irY, \(N + l,k,{p} n ^\H T \±,{p} n )\ 2 

n,pj ,ki 

x 8{E k - E p - u±)p(j3H ), (11) 



r 2± = 2n J2 \{N + l,p,{k} 2n -i\H T \±,{kh 



x 5(E p -E k -u ± )p{PH ), 



(12) 



where T± is decay rate for the excited |+) or ground state 
|— ) of the qubit, and u± is defined in Eq. (J2J. State 
| + ,{p}„), for example, denotes the excited state of the 
qubit with n quasiparticles in the reservoir with energies 



E p = J el + A2 



\+,{p}n) = |+>®|P1>- 



,Pr, 



(13) 



The state |+,{fc}2n) denotes the excited state of the 
qubit with n broken Cooper pairs in the box, lead- 
ing to the appearance of 2n quasiparticles with energies 



AT 



\ + ,{k}2n) = 1+,**,. 



(14) 



In the following, we concentrate on the decay rate of 
the qubit excited state in the open system, i.e., evalu- 
ate r^ p . In order to calculate this decay rate we take 
into account one-electron processes in the lowest order in 
quasiparticle density. We assume that in the first con- 
tribution to T°F, Eq. Hll|) . all quasiparticles are in the 
reservoir and one of them is tunneling into an unoccupied 
state of the CPB; in the second contribution, Eq. (|12f> . all 
quasiparticles are in the box and one of them is tunnel- 
ing out into an unoccupied state of the reservoir. Keep- 
ing this in mind, the density matrix for first process can 
be reduced by tracing out irrelevant degrees of freedom: 
p op (f3H ) — Tr^p(PH ), and Eq. (fTT|> becomes 



r° p = ttJ2 \{N + l,k,{p} n -i\H T \N,ip} f 



x 8(E k - E p - uj±)p op ((3H ). 



(15) 



Taking into account that only one quasiparticle is trans- 
ferred through the junction by the action of Hamiltonian 
Ht, and performing the sum over momenta in Eq. 115fl . 
the contribution of the first mechanism is simplified to 

T? P = t£ I 1, M H T \N,pi) \ 2 6(E k - E Pl - w+) 



x exp 



E 



Pi 
T ) 1 



(16) 



where the exponential factor is the low temperature 
(T <C A r ) approximation of the Fermi function. The ma- 
trix element (N + 1, k\ Ht \N,p) can be calculated using 
the particle conserving Bogoliubov transformation^ and 
is equal to (N+ l,k\H T \N,p) = 2(t pk u p u k - t kp v p v k ), 
where u p , v p are coherence factors 



<=2 



E, 



and v 2 p = - 



1 - 



E, 



By changing the sum to an integral in Eq. (|l(jfl and inte- 
grating over E k , we get the following expression for r° p : 
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r op = 9_ f°° dE <d(E p +u; + -A b ) 

x (E p (E p +oj + )-A r A b )exp(~^-y (17) 

where 0(x) is the step function, g = is dimensionless 
conductance. R is resistance of the tunnel junction in the 
normal state 

R- 1 =4ne 2 J2\t P k\ 2 S(e P )S(e k ). 

p.k 

. i 



Assuming that mismatch between superconducting gap 
energies in the box and reservoir is small, A r — A b +cu + > 
0, which corresponds to most charge qubit experiments, 
expressions for T° p can be simplified. The leading con- 
tribution to the decay rate at low temperatures given by 
the first mechanism is equal to T° p = W(lo+ 1 A r , A b ) 



W(lo + , A r ,A b ) 



gy/A~; {A r -A b +u}+) 
8\/27r v / A r + A b +uj + 



exp 



A r 



-u 



r.2,- 



T 



, (18) 



r 



where U (a, b, z) is the confluent hypergeometric function. 
At low temperature (T <C A r — A b + lu + ) asymptotic 
result for W(u>+, A r , A b ) is simply given by 



W(LO+ ,Ar,A b 



gVK/F A r -A b 



4V2tt 



■ exp 



A r 



As expected, the decay rate due to the first mechanism 
is exponentially suppressed due to the fact that it costs 
energy A to bring a quasiparticle from the normal parts. 

The contribution of the second mechanism given by 
Eq. (|12[) depends on the density matrix of the box. The 
initial state of the qubit corresponds to the even-charge 
state in the CPB. Statistical weight of the states with 
even number of quasiparticles in the dot, 2, 4, 6, 2n, is 
determined by the density matrix p2n(PHo)), 



exp 



P2n(/3H )) = Tr {p}P (l3H ) 




(2n)\Z c 



(20) 



where Z cv — cosh (zb(T, 5 b )) is the partition function for 
the dot with even number of electrons^, and z b (T, S b ) is 



Zb(T, 5 b ) 



E ex p(-|^ 

k V 



T A b 
A fc b b 



Ab 
T 



(21) 

According to Eq. I|12|l and Eq. I|2()|l. the contribution to 
the decay rate due to the second mechanism is obtained 
by summing over the states with even number of quasi- 
particles with appropriate statistical weight 

T° P = 27T Y, l(^ + l,P.W2«-l|ffr|+,{fc}2n>| a 

n,p,kj 

x 6(E p - E kl - cj+)2np 2n {[3H ), (22) 



where, for example, (N + l,p, {fc}2n— 1| is a state cor- 
responding to the charge on the box equal to le, 2n — 1 
quasiparticles in the box, and 1 quasiparticle in the reser- 
voir 

\N + l,p, {k} 2n -i) =\N + l,kx,...,kj,..., k 2n -i)®\p) . 

(23) 

The additional factor of 2n in Eq. I|22l) is the result of 
the summation of 2n identical terms in Eq. (|12|) . The 
tunneling matrix element in Eq. I|22l) is determined using 
the particle-conserving Bogoliubov transformation and is 
dependent only on p and k\. Therefore, by doing the sum 
over the other momenta kj, one gets the following result: 

17 = tt^KW +l,p\H T \N + 2,k 1 )\ 2 S(E p - E kl -Lu + ) 

p,ki 



x exp 



T 



[z b (T,6 b )] 2n - 
(2n-l)!Z 



(24) 



Changing the sum to an integral and integrating over E p , 
we obtained the following expression: 



(E k (E k + uj+)- A r A b ) 



^m-AD 



E k 



x e(E fc +w+-A r )exp l-^jtanh(2 6 (T,<5 6 )) .(25) 

The integration can be performed assuming that mis- 
match between superconducting gap energies in the box 
and reservoir is small, and A b — A r +uj + > 0. Comparing 
Eqs. I|17l) and l|25|l . one notices that the answer for r°, p 
can be expressed via by permuting A r <-> A b 

r° p ( w+ , A r , A b ) = r° p (w+,A b ,A r )tanh (z b (T, S b )) . (26) 

Taking into account Eq. (|18|l and Eq. (|2(j|l , we find the 
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total quasiparticle decay rate for the open system 

r,£ +1)H+) = W(uj+,A r ,A b ) + W(w + , A b ,A r ) 

x t&nh {z b (T, S b )) , (27) 

where the first term corresponds to the first mechanism 
and is dominating for the open system. The simplified 
are discussed in Sec. IV. 



results for 



-1+) 



Decay rate evaluation for the ground state of the qubit 
|— ) can be done similarly provided that the final state 
\N + 1) of the qubit is lower in energy than the initial 
state I—). 



III. STATES OF THE QUBIT IN AN ISOLATED 
SYSTEM 



A. Thermodynamic properties of the qubit with 
fixed number of electrons 

When the number of electrons in the qubit is fixed, 
parity effects become important at low temperatures, 

T < T*, T b *, -■ ■ - ■ - 

equal to 



where T b * is defined in Eq. J3J) and T r is 



Tr ln(A r /<5 r ) ' (28) 

Here, 8 r is the mean level spacing in the reservoir. The 
density of quasiparticles in the qubit with even number 
of electrons is small, n oc exp (— 2A/T), since at low tem- 
peratures all electrons arc paired, and it costs energy 2A 
to break a Cooper pair. In the odd-charge state, an un- 
paired electron is present in the system even at zero tem- 
perature. It is important to understand the probability 
of finding a quasiparticle in the CPB since the qubit will 
not work if an extra electron is present in the box. To 
find out whether it is favorable or not for a quasiparticle 
to reside in the CPB, we calculated the difference in free 
energy AF = Fi—Fq between two states: with and with- 
out a quasiparticle in the box (iq and Fq, respectively). 
At the operating point, the free energy difference for the 
qubit with even (AF ev ) and odd (Ai^dd) total number 
of electrons is given by the following expressions: 



AF nv — —E a 



^-Tln[tanh(z b (TA))] 



Tln[tanh(z r (T, 6 r ))} 



AFr 



(1(1 



Ej_ 
2 

Tln[coth(z r (T, 5 r ))] 



Tln[tanh(z b (TA))] 



(29) 



(30) 



Negative value of AF indicates that free energy is lower 
for a quasiparticle in the CPB. Using these expressions, 
we can calculate thermodynamic probability P(T) to find 
an unpaired electron in the box as a function of the phys- 
ical parameters 



P(T) = 



1 



Z 



i + Z Q exp(AF) + l' 



(31) 




T(K) 



FIG. 2: (color online). Main panel: temperature dependence 
of the number of quasiparticles in the box at the operating 
point (N g = 1). Dash-dot (black) line corresponds to even 
number of electrons, solid (blue) line - odd number of elec- 
trons, dash (red) line - open system. Physical parameters 
are chosen in correspondence to typical qubit experiments: 
A r = A b = 2.4K, E c = 0.25K, Ej = 0.3K, T 6 * = 210mK, 
and T* = 160mK (see Eqs. and for definition of T 6 * 
and T*). Inset: temperature dependence of the number of 
quasiparticles in the CPB in the vicinity of T*. 



where ^i(o) = ex P ( — /3-F\(o)) is the partition function 
with one (zero) unpaired electrons in the box. The ex- 
pression for the free energy difference AF op for the open 
system (fixed chemical potential regime) can be obtained 
using A_F CV and taking the limit of infinite volume of 
the reservoir (S r — » 0). Temperature dependence of the 
quasiparticle probability in the CPB for realistic param- 
eters is plotted in Fig. 2. 

As shown in Fig. 2, at high temperatures T > T b * the 
probability of having an extra electron in the CPB coin- 
cides for an open and isolated qubit. At this tempera- 
ture the number of thermal quasiparticles in the system 
is large, and parity effects are not important. Parity ef- 
fects start to manifest themselves below the characteristic 
temperature T fc * , when the number of quasiparticles in the 
box is of the order of unity. As can be seen from Fig. 2, at 
the temperature T* the probability of having a quasipar- 
ticle in the CPB is negligible in the even-charge state, as 
well as in the open system. In the case of odd charge state 
of a qubit, lowering the temperature enhances quasiparti- 
cle poisoningi^. This effect can be explained as a compe- 
tition of two contributions to the free energy: a charging 
energy gained by tunneling to the box and entropy con- 
tribution proportional to the ratio of the volumes of the 
reservoir and box, ~ V r /V b . At the temperature T s 



T, 



E c - Ej/2 



In 



(ft) 



(32) 



the entropy contribution becomes smaller, and the quasi- 
particle resides in the CPB. Thus, for the odd-charge 
state there is only a certain intermediate temperature 
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range when qubit can work, i.e. can be prepared in the 
"good" quantum state. For physical parameters used in 
Fig. 2 T s is approximately equal to 20mK. 

The presence of a quasiparticle in the box can be 
studied experimentally by measuring the periodicity of 
Coulomb staircase According to the ob- 

tained results, the Coulomb staircase for an open system 
should be 2e periodic below the temperature T b . The 
qubit with fixed number of electrons should have two 
distinct types of behavior corresponding to even and odd 
total number of electrons in the box and reservoir. In the 
former case the Coulomb staircase is similar to the one in 
the open system, while in the latter Coulomb staircase is 
le periodic for temperatures above T b * , then 2e periodic 
from T fc * to T s , and again le periodic for T < T s . 

In principle the probability to find a quasiparticle in 
the CPB can be lowered, and the qubit can be brought 
into the desired quantum state. We discuss several ways 
of doing this in Sec. IV. However, even if the quasiparti- 
cle is in the reservoir at the initial moment of time, once 
the qubit is excited and quantum oscillations emerge, the 
time of the oscillations is determined by the kinetics, i.e., 
by quasiparticle tunneling rate. 



B. Quasiparticle decay rate in the isolated system 

Let us turn to the discussion of the lifetime of the 
charge qubit in the regime with the fixed number of elec- 
trons. In order to calculate the quasiparticle decay rate, 
we proceed in the same manner as in the open system. 
The decay rate for the even number of electrons is cal- 
culated by averaging over initial states with even parity 
density matrix p2 n {^H( S ). This situation corresponds to 
having an even number of electrons in the box and reser- 
voir. The appearance of quasiparticles in the system oc- 
curs at the expense of breaking Cooper pairs. Using the 
results of analogous calculation in Eq. I|2t)l) , we can write 
the expression for the total decay rate 

r |^+i)H+) = W>+,A r ,A fc )tanh(z r (TA)) 

+ IU(w+,A b ,A r )tanh(z 6 (T,<5 b )), (33) 

where W(u + , A r , A b ) is defined in Eq. (jT%|l . The first 
term here corresponds to the first mechanism given by 
Eq. 111|) and averaged over the even-parity initial state. 

In the odd-charge case the decoherence rate is the 
largest since a quasiparticle is present in the system at 
T = 0. Initial configuration of the system corresponds 
to having an odd number of quasiparticles in the reser- 
voir. The reduced density matrix for this initial state 
p 2n -i(/3H ) is then given by 



exp 



P2n - 1 (/3H ) = Tr {k} p(/3H ) 




where Z dd = sinh (z r (T, 6 r )). Using Eq. I|ll|) . we write 
the contribution to the decay rate of the first mechanism 

rf d - 2tt ]T |(^ + i,Mrf 2 „- 2 |tf T |+,M2„-i)| 2 

n,pj ,fc 

x 5{E k - E P1 - uo+){2n - l)p 2n _ 1 ((3H ). (35) 

Going through the same arguments as in Eq. I|24|l . the 
expression for r° dd can be simplified 

IT" = kY,\(N + l,k\H T \N, Pl )\ 2 8{E k -E pl -u J+ ) 

p lt k 



. ( , M) ,_^WWrA)] 



2n-2 



T J ^ (2n-2)!Z odd 

/ n— 1 v ' 

Summing over Ek, we get 



(36) 



podd _ _9_ 

1 4tt 



dEL 



(E P (E P + lo+) -A r A b ) 
{{E p +u + y-Al){El-AD 



E, 



(2n-l)!Z odd 



(34) 



x e(£ p +c + -A b )cxp(-^j coth (z r (T,5 r ))(37) 

Taking into account Eq. JT7|) and Eq. JTB), T1 dd is equal 
to 

T° dd = W(uj+, A r , A 6 ) coth (Zr(T,5 r )) ■ (38) 

In order to find the contribution of the second mecha- 
nism, one has to average over initial states of the CPB. 
The initial configuration of the box corresponds to the 
even-charge state and is the same for open and isolated 
qubits. Therefore, contribution of the second mechanism, 
r^ dd , is given by Eq. J2EJ- 

Total decay rate Fj^^^^ with odd number of elec- 
trons in the qubit is the sum of r? dd and T^ dd 

r|w d +i)H+} = ^(^ + ,A J .,A b )coth(z J .(T, ( 5 r )) 

+ IU(^ + ,A fc ,A r )tanh(2 b (r,,5 fc )) (39) 

At low temperatures, T < T r *,T fc *, the first term here 
is dominant as z(T, S) 1. The detailed analysis of 
the low temperature asymtotics for different experimen- 
tal regimes is presented below. 



IV. DISCUSSION OF THE RESULTS 

Temperature dependence of the quasiparticle decay 
rate for different experimental realizations of the qubit 
is shown in Fig. 3. As it is clear from the figure, at exper- 
imentally relevant temperatures T «T r * , T b * , the largest 
decay rate corresponds to the odd-charge case. In the 
vicinity of T* defined in Eq. J5SJ), the decay rate is grow- 
ing quickly due to the appearance of a large number of 
quasiparticles in the reservoir. As we approach the tem- 
perature T — T fc * given in Eq. @, which corresponds to 
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FIG. 3: (color online) Temperature dependence of the quasi- 
particle decay rate. Dash (red) line corresponds to the open 
system, dash-dot (black) - even number of electrons, solid 
(blue) - odd number of electrons. Here, we used the same 
physical parameters as specified in Fig. 2. 

the appearance of quasiparticles in the Cooper pair box, 
parity effects become irrelevant and decoherence rates for 
different cases coincide. 

Results obtained in Sec. II and III allow us to quantita- 
tively estimate the decoherence rate due to the presence 
of quasiparticles in the system. For simplicity we assume 
that superconducting gap energies are the same in the 
box and reservoir, A& = A. r = A, and temperature is 
low, T < T*, T£ <C to+ , corresponding to typical qubit 
experiments (cj + is defined in Eq. ©)■ In this approxi- 
mation, for an "open" qubit the decay rate is 



r op av j " / ^+ I I 

i |7v+i)^|+> - 8 0F V A + ^+/2 V T) 

+ ^VAT^r xp (-^)- (40) 

The leading contribution to the decay rate for the open 
system r°^ +1 ^| + ^ is proportional to exp(— A/T) since 
energy equal to A is required to bring an electron from 
the normal parts. 

In the even-charge case the leading contribution to 
rf7v + i)^| + ) is exponentially small, oc exp(— 2A/T), in 
accordance with the energy necessary to break a Cooper 
pair, 

gT I u; + A / 2A\ 
V+DH+) " A+c + /2 6V P \~TJ 

80FV A+cu+/2 5 b 1 \ T ) ' 

Note that Eq. I|41|) and the last term in Eq. (|4U[) are essen- 
tially the upper bounds for the contributions to the decay 
rate coming from unpaired electrons which originate in 
the isolated qubit and in the Cooper pair box, respec- 
tively. Indeed, we assumed in the derivation of Eqs. 14U|) 



and (|41|) that the decay rate is limited by the thermo- 
dynamic probability of the unpaired state, without dis- 
cussing the kinetics of pair-breaking leading to such state. 
It is clear from Fig. 3, however, that the above mentioned 
contributions are negligibly small at low temperatures. 

The decay rate of the qubit with an odd number of 
electrons is much larger 

F odd _ jrfr_ I ^+ 




The leading contribution to the decay rate is temperature 
independent since the number of quasiparticles is finite 
even at T — 0. According to the Eq. I|42|l . the lifetime of 
qubit in the odd-charge case is determined by the conduc- 
tance of the tunnel junctions and mean level spacing of 
the reservoir. For typical experiments quantum conduc- 
tance, g is of the order of one; S r depends on the volume 
of the reservoir and varies between 10 -10 and 10~~ 12 eV. 
With these parameters, decay rate rj 5 ^ can be es- 

timated as 10 5 — 10 3 Hz consequently. This is a substan- 
tial contribution to the decoherence of the isolated charge 
qubits, which limits qubit operation on a fundamental 
level. However, this decay rate is much smaller than the 
present estimates for decoherence in charge qubits; see, 
for example, the recent review by Devoret et. ai&. 

We assumed so far that in the initial state an unpaired 
electron resides in the reservoir and finally (after the re- 
laxation) ends up in the box. Then, one can tune the 
qubit to the charge degeneracy between le and 3e and 
still have charge oscillations for some time until the quasi- 
particle escapes into the reservoir. In this case, the quasi- 
particle escape rate can be calculated in a similar way and 
is proportional to the conductance of the tunnel junction 
and level spacing in the Cooper-pair box: r odd cx g5b- 

In principle, quasiparticle poisoning can be decreased 
by tuning the proper parameters of the system such as 
superconducting gap energies A r! b, charging and Joseph- 
son energies, temperature, volumes of the box and reser- 
voir. For example, it can be done by adjusting gap ener- 
gies A r b with the help of oxygen doping^ or magnetic 
fieldii. The latter is easy to implement since magnetic 
field is already used in charge qubits to tune Josephson 
energy Ej. In addition to the suppression of gap ener- 
gies, large magnetic field (H > H c i) can create vortices 
in the reservoir. The vortex acts as a quasiparticle trap 
since an unpaired electron gains gap energy A by residing 
inside the vortex. 



V. CONCLUSIONS 

We demonstrated that the presence of quasiparticles in 
the superconducting charge qubit leads to the decay of 
quantum oscillations. Two experimental realizations of 
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charge qubit are considered here, corresponding to open 
and isolated system (in the former, the number of elec- 
trons is not fixed). Once the qubit is excited and quan- 
tum oscillations emerge, the decay of these oscillations 
is determined by the quasiparticlc tunneling rate to the 
Cooper pair box. We calculated temperature dependence 
of the quasiparticle decay rate in the charge qubit. This 
decay rate is exponentially suppressed in the open sys- 
tem as well as in the isolated system with an even number 
of electrons. However, in the case with an odd number 
of electrons in the system, the quasiparticle decay rate 
is not exponentially suppressed and is estimated to be 



10 5 — 10 3 Hz depending on the volume of the supercon- 
ducting reservoir and conductance of the tunnel junc- 
tions. 
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